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! , ABSTRACT: The Hadron Resonance Gas (HRG) model is considered to study the QCD equation of state 

for the case of nonzero external magnetic fields. Thermodynamic observables including the pressure, energy 
i density, entropy density, magnetization and the speed of sound are presented as functions of the temperature 
and the magnetic field. The magnetization is determined to be positive, indicating that the hadronic phase 
i of QCD is paramagnetic. The behavior of the speed of sound suggests that the deconfinement transition 
temperature is lowered as the magnetic field grows. Moreover, a simple correspondence is derived, which 
^ relates the magnetic catalysis of the quark condensate to the positivity of the /3-function in scalar QED. 
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1 Introduction 

Quantum Chromo dynamics (QCD) is the established theory of the strong interactions. QCD predicts that 
strongly interacting matter exhibits at least two very different forms: the hadronic phase at low temperatures 
T and the quark-gluon plasma (QGP) phase at high T. Whereas the relevant degrees of freedom are color- 
singlet objects in the hadronic regime, they are colored particles in the quark-gluon plasma state. The 
transition between the QGP and the hadronic phase took place as the early universe expanded and cooled, 
about 10~ 5 seconds after the Big Bang. On the other hand, the same transition is also reproduced in heavy- 
ion collisions conducted in contemporary experiments at e.g. the Relativistic Heavy Ion Collider (RHIC) and 
the Large Hadron Collider (LHC). A particularly interesting aspect of the transition is how thermodynamic 
observables change as the system passes through the region separating the hadronic and QGP phase. The 
relations between these observables constitute the equation of state (EoS) of the system in equilibrium. 

A striking observation made in heavy-ion collision experiments is that the flow of strongly interacting 
matter can be adequately described in terms of nearly-ideal relativistic hydrodynamic models, see, e.g. 
Refs. [1-4]. The dependence of the EoS on state parameters like the temperature T and the chemical 
potential [x are necessary input to these models. For a non-central heavy-ion collision, however, an external 
magnetic field B is also generated by the spectators. Since the strength of this magnetic field reaches up to 
the hadronic scales [5-7], it can have a significant impact on the properties of the transition and of the EoS. 
Similarly strong magnetic fields are expected to be present in dense neutron stars [8] and to be generated 
during the electroweak transition in the early universe [9]. In each of these cases, the interplay between the 
strong dynamics and the coupling to the external field can induce new and exciting phenomena. Examples 
include the the chiral magnetic effect [10, 11] or the decrease of the transition temperature with growing 
magnetic field [12, 13]. 
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Thus, a clear theoretical understanding of the dependence of the EoS on T, \x and B is desired. Most 
of our knowledge about the B = \i = EoS comes from lattice Monte-Carlo studies, see e.g. Refs. [14- 
16]. The inclusion of a finite chemical potential \i poses conceptual problems, when it comes to the lattice 
approach. On the contrary, nonzero magnetic fields are straightforwardly simulated using standard Monte- 
Carlo algorithms. Still, lattice results about the EoS for nonzero external magnetic fields are yet absent 
from the literature. 

On the other hand, there is another, remarkably simple approach, which can be used to access the low- 
temperature region of the EoS: the Hadron Resonance Gas (HRG) model. Within this model, the hadronic 
phase of QCD can be studied, even at nonzero chemical potentials or external magnetic fields. Comparison 
to lattice QCD results reveals that the HRG description gives a good approximation of thermodynamic 
observables even up to temperatures just below the transition region, both at zero chemical potential (see, 
e.g., [16-18]) and at nonzero chemical potential (see, e.g., [19-22]). Still, the HRG description has not yet 
been employed for the case of nonzero magnetic fields. 

In this paper, the QCD EoS for nonzero magnetic fields is determined within the HRG model. Sec. 2 
is devoted to the discussion of thermodynamic relations in the presence of a magnetic field. In Sec. 3, the 
HRG approximation is detailed, and the free energy for individual hadrons is calculated. The contributions 
from each hadron are summed up in Sec. 4 to obtain the total free energy density, from which the whole 
equation of state is reconstructed. The B- and T-dependence of thermodynamic observables including the 
pressure, energy and entropy density, magnetization and speed of sound are calculated. Finally, in Sec. 5, 
we conclude. 



2 Thermodynamics in an external magnetic field 

The quantity on the top of the hierarchy of thermodynamic relations is the thermodynamic potential - 
which we refer to as free energy. In terms of the partition function of the system, this free energy reads 
T = —TlogZ, and in the presence of a constant, external magnetic field B is written as [23, 24] 

F = £ - TS - BMb, (2.1) 

with £ the energy 1 , S the entropy and Mb the magnetization. These observables satisfy the differential 
relations 

8T=- S > dB=~ MB > W = - P - (2 ' 2) 

We also define here the corresponding densities as e = £/V, s = S/V, f = F/V and = Mb/V. In the 
thermodynamic limit, V — >• oo, differentiation with respect to the volume simplifies to a multiplication by 
1/V, and thus the pressure is given by 

P = -^ = -/ = -(/ vac + / therm ), (2.3) 

where we anticipated that the free energy separates into a vacuum and a thermal contribution, see Eq. (3.6) 
below. From Eqs. (2.1) and (2.3) the energy density is calculated as 

e = Ts + Bm B -p. (2.4) 

Another observable of interest is the speed of sound c s , which is defined in terms of differentials at constant 
B, 
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(2.5) 



1 Note that E in Eq. (2.1) denotes the total energy of the system, with the work necessary to maintain the constant external 
field also taken into account [23]. 
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We remark that the definition (2.3) leads to a pressure that is isotropic in space - solely due to the assumption 
that the free energy J- is an extensive quantity, and (in the thermodynamic limit) is proportional to V . As 
a result, the speed of sound will also be isotropic. Defining p in terms of the diagonal elements of the stress 
energy tensor leads to anisotropic pressures (see, e.g., Ref. [25]). However, this anisotropy is canceled, if the 
work done against the Lorentz force due to the magnetization is also taken into account [26], resulting in a 
pressure that coincides with the thermodynamic definition of Eq. (2.3). 

3 Free energy density in the HRG model 

In the HRG model [27], the thermodynamic potential, Eq. (2.1), of the system is approximated in the 
thermodynamic limit, V — > oo, by the partition function of a gas of non-interacting free hadrons and 
resonances [28]. Thus, the free energy density of the model is written as the sum of independent contributions 
coming from non-interacting hadrons h, 

f = ^2 d h- fh ({eB,T},{m h ,q h /e,s h ,g h }) , (3.1) 
h 

where each contribution fh depends on the external parameters (the magnetic field in elementary charge 
units eB and the temperature T), and internal properties of the hadron (mass rrih, spin Sh, charge qh/e 
and gyromagnetic ratio gh)- Each hadron enters the sum with a certain multiplicity d^. The hadrons taken 
into account extend from pions up to the S° baryon, as listed in the latest edition of the Particle Data 
Book [29], and are tabulated in Table 1. The experimental values for the gyromagnetic ratios in Ref. [29] 
are known only for a few hadrons, and only with large uncertainties (except for the proton and the neutron). 
Therefore we decided to take the gyromagnetic ratios to be gh = 2qh/e, as dictated by universal tree-level 
arguments [30]. This corresponds to the assumption that the considered hadrons are point-like objects; e.g. 
neutral hadrons have gh = 0. A possible improvement of the method is to take into account the correct 
gyromagnetic ratios. The energy levels for s = 1/2 and s = 1 particles with anomalous magnetic moments 
have been discussed in Refs. [31, 32]. These, however, lead to more complicated expressions for the free 
energies, for which the Landau sums cannot be performed in general. 
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Table 1. List of hadrons and resonances taken into account in the HRG description. 



3.1 Energy levels 

To reconstruct the free energy in the low-temperature region using Eq. (3.1), contributions from each particle 
type are summed up, with the assumption that the interaction between them is negligible. Let us therefore 
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take a free relativistic particle with momentum p = (p x ,Py,Pz) aiK1 mass m, in the presence of a magnetic 
field of magnitude B, pointing in the positive z direction. We consider the particle to have spin s and charge 
absolute value q (thus, in our notations qB is always positive). The component of the spin in the direction 
of the magnetic field is a conserved quantity, and can assume the values s z = —s, —s + 1, . . . s. The Landau 
levels are labeled by the index k. With these notations, the energy levels of a charged particle {q > 0) in 
the presence of the magnetic field are given as [33], 

E(p z , k, s z ) = ^pl+m 2 + 2qB(k + l/2-s z ), (3.2) 

while the energy levels for the neutral particle {q = 0) are 

£ (p) = Vp 2 + m 2 . (3.3) 

As mentioned above, a generalization for the description of anomalous magnetic moments for s = 1/2 and 
s = 1 has been developed in Refs. [31, 32]. In that approach, using a magnetic field-dependent gyromagnetic 
ratio, certain inconsistencies in the s = 1 theory can be resolved. On the other hand, the description of the 
spin-3/2 theory in terms of Rarita-Schwinger fields [34] is known to exhibit non-causal behavior [35, 36]. 
Whether the simple formula (3.2) adequately describes the dispersion relation of a spin-3/2 particle is 
therefore not obvious. In fact, we will show in Sec. 3.5 that - unlike the other spin channels - the s = 
3/2 sector would give a negative contribution to the pressure for any nonzero magnetic field, marking an 
instability in the theory. For this reason we do not consider resonances with s = 3/2 or higher in the model. 

At arbitrary finite temperature, the free energy density (at vanishing chemical potentials) for a charged 
particle can be written as (see Ref. [37] for the B = relation, and e.g. Ref. [38] for the analogous expression 
at B ^ 0), 

s z k=0 J ^ ' 

where the lower sign corresponds to bosons (s integer) and the upper one to fermions [s half-integer) 2 . The 
same for a neutral particle, on the other hand, is given by 

/nM = =F£ /Kf + ^li^)). (3.5) 

We separate the total free energy density into a vacuum (T = 0) and a thermal part as 

r c (s) = f( S )\ T=Q , / th - m ( s ) = f( s ) - r c ( S ). (3.6) 

The vacuum terms in Eqs. (3.4) and (3.5) are ultraviolet divergent and need to be regularized. We use 
dimensional regularization with d = 1 — e. After separating the divergent contribution, the renormalization 
of the free energy density is carried out by subtracting the B = term, and by performing the renormalization 
of the pure magnetic energy B 2 /2. As we will see, the latter is equivalent to renormalizing the elementary 
electric charge e [39-42]. The vacuum free energy density has been obtained previously in e.g., Refs. [38-49], 
however, the connection to electric charge renormalization has not been stated in all cases. In fact, it is a 
remarkable feature of the background field method, that electric charge renormalization - and, accordingly, 
the coefficients of the /3-function - can be determined solely from how the free energy depends on an external 
magnetic field [41]. In particular, this approach also gives insight into how the /3-function can be related to 
the mass-dependence of the free energy - and, thus, to the magnetic catalysis of the quark condensate (see 

2 Antiparticles are not taken into account here. 
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Subsec. 3.4). We remark that our renormalization prescription ensures, that the contribution to the free 
energy density diminishes as the mass of the particle increases, and thus the sum of Eq. (3.1) in each spin 
sector will be convergent. 

On the other hand, the thermal part is explicitly finite (reflecting the fact that each divergence is 
independent of the temperature), and can be determined by numerical integration and summation. 

3.2 Vanishing magnetic field 

In order to carry out the renormalization of the B > free energy density, first it is necessary to determine 
the B = contribution. The vacuum free energy density ati? = 0in<i = 3 — e dimensions is given by 

/-c (S)jB = ) = T I(2 s + l).^y _^-!p_7^T^, (3.7) 

where the lower sign is for bosons, and the upper for fermions. Here the scale ji appeared to fix the dimension 
of the above expression to four. To relate this integral to the nonzero B case - following Ref. [46] - we 
rescale the momenta as p — > p^2qB, with qB being an arbitrary dimensional scale. We denote x = m? /2qB 
and perform the integration using the formula (A.l) to obtain 

f™(s, B = Q) = T 1 (2s + 1){2qB) 2_± ^ r( _ 2 + e/2 ) X W2_ (3. 8) 



We expand in e using Eq. (A. 2), 



H log k log(x) 

e 4 2 2 & \4irfi 2 J 2 BK J 



(3.9) 



8vr 2 

Expressing this with x = m 2 /2qB cancels all ^-dependence, of course. 
3.3 Nonzero magnetic field 

Let us now consider a particle with charge q and spin s in a magnetic field B. The energy levels are given by 
Eq. (3.2). In terms of these levels, the vacuum free energy density is written using dimensional regularization 
as 

/ vac («) = T^EE|f-^ / ^fe Vfl + ™ 2 + 2*B(k + 1/2 - s z ), (3.10) 
k=o s z J { ) 

where again the upper sign corresponds to fermions and the lower to bosons, and \x is the scale related to 
dimensional regularization. We abbreviate a = 1/2 — s z and integrate in p z (using formula (A.l)) to get 

r c {s) = ± { ^{^) e/2 r(-l + e /2)^C(-l + e /2,* + a), (3.11) 

where the sum over k was converted into a Hurwitz-^ function, see Eq. (A. 3). Now expanding in powers of 
e, using Eq. (A. 2), we obtain 

/-« ? {{-h * - - >) H - ^ + - + < 

(3.12) 

To calculate the change in the free energy density due to the magnetic field, we subtract the B = 
contribution, Eq. (3.9), 

AT-M - ^ E [(! - t - * . {^) + (i ~\ + f ) " f(- W •) " M *M 

(3.13) 
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where we used that Y^a( a ~ V^) = 0- 

This expression is still divergent, as it contains the purely magnetic field-dependent term ~ B 2 /e. In 
order to cancel this divergence, we have to redefine the free energy density by including in it the energy 
density of the magnetic field, B 2 /2 [39-42]. The divergence is then absorbed into the renormalization of the 
electric charge, and simultaneously, into that of B, 

A/ vac, = A/ vac + ^ ) B 2 = Z g B 2 r , q 2 = Z^q 2 , QrB r = qB , (3.14) 

with the renormalization constant 

^Tj|£0-,->o g (^))(i- a + 4 (3.5) 

which, for a spin-1/2 particle, 

Zr°* = 1 + ^i Pin ° r (" \ + 1 + ^ ) > Z 3 ™ = ^ (3-16) 

reproduces the well-known expression 3 (see, e.g., Ref. [50]), with the leading coefficient of the spinor QED 
/3-function appearing in front of the divergence. Similarly, for a spin-0 particle we get 

^ calar = 1 + 2^ q r (~~ + 7 + lo § {0^) ) > ^ = ^2 » (3-17) 

where the scalar QED /3-function coefficient enters (cf. Ref. [41]). Here m+ = m is a constant, which is fixed 
to the physical mass of the particle (see discussion in Subsec. 3.4). The renormalization of Eq. (3.14) then 
leads to 



A/ vac,r (s) = _r ^ J- x + a ) _ _ l og(x ) + + (l og(x ) + 1 



(3.18) 



which can be rewritten in a compact form in the original variables, 



A/ vac ' r (s) = — =f 



B 2 ^(qB) 2 _ . 



2 8?r 



2 



C(-l,s + l/2-a,) 



+ i 2, + i ) ■ ( t " T log(x) + lQg( 2zl + 1 0- - M« + 1)) 



(3.19) 



We remark that instead of using dimensional regularization, the free energy density in the s = and 
s = 1/2 sectors can also be calculated using Schwinger's proper time formalism [39], which gives the same 
results [41, 44]. 

Note that the renormalization constant Z q is such, that the renormalized A/ vac ' r at 0(B 2 ) is exactly 
B 2 /2 - in accordance with the expectation that the free energy, at T = and at small magnetic fields, comes 
exclusively from the external field itself. Indeed, using the asymptotic behavior, Eq. (A. 5), of the Hurwitz £ 
function, it is easy to check that the only contribution to A/ vac,r at 0{B 2 ) is the pure magnetic term, and the 
small B (large x) limit of the square parentheses in the second contribution of Eq. (3.18) is zero. Note that 
this also ensures, that the large mass limit of the matter contribution vanishes. Different renormalization 
schemes have also been used in the literature. For example, the scheme employed in Ref. [46] corresponds 

3 The factor 1/2 in front of /3^" no1 [ s canceled, if the contribution of the antiparticle is also taken into account. To convert 
to cutoff regularization, see Eq. (A. 7). 
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to dropping all mass- independent terms in Eq. (3.13), whereas in the one used in Ref. [48], an additional 
~ {qB) 2 Q'{— 1, a) finite term is also subtracted. These schemes are connected to that of Eq. (3.15) by a 
finite renormalization, which would produce a renormalized vacuum free energy of the form B 2 /2 + 0(B 2 ) 
- instead of defining the total quadratic term to be B 2 ./2. Accordingly, the different renormalization in 
Refs. [46, 48] leads to a A/ vac,r (and thus, also to a magnetization mg) that grows logarithmically with the 
mass - instead of approaching zero for m — > oo, as expected on physical grounds. The prescription (3.15) is 
the only choice, for which the total quadratic term is B 2 /2, and the large mass limit of the free energy (and 
the magnetization) vanishes. 

If more particles (possibly with different masses) are present in the system, the renormalization constant 
Z q is extended to absorb the divergences coming from the interaction of each particle with B. The pure 
magnetic energy B 2 /2 will be unchanged, independently of the number and properties of these particles. 

3.4 Renormalization and magnetic catalysis 

Let us consider the mass-dependence of the vacuum free energy density. In this respect, one has to care- 
fully distinguish between the actual mass m of the hadron, and the fixed mass m* that appears in the 
renormalization prescription, Eq. (3.15). Writing out the variables explicitly, the first relation of Eq. (3.14) 
reads 

Af ac < r (m, m*) = A/ vac (m) + ^ Z,(m*). (3.20) 

The renormalization ensures that the total renormalized free energy density up to 0(B 2 ) equals the pure 
magnetic contribution B 2 /2 for a hadron of mass m = m+, 

Ar^(m,m,)\ n = ^ + 0(B% (3.21) 

The renormalization prescription (i.e., m*) does not change, if the mass of the hadron is varied - in this 
sense the employed scheme is mass-independent. 

Let us now consider the quark condensate at T = 0, which is given in terms of the vacuum free energy 
density as 

^• s --kr- (3 - 22) 

where m q is the mass of the quark. We can write the change in the condensate due to the magnetic field 
through the hadron sigma terms (d(m 2 ) / dm q ) as 



9Ar- r (m,m,) d{m 2 ) = 9A/ vac (m) d(m 2 ) 
d(m 2 ) drriq d{m 2 ) dm q 



where we inserted Eq. (3.20). One can easily check using Eq. (3.13) for A/ vac , that the result is a positive 
condensate Aipijj > of 0(B 2 ), in agreement with the well-known magnetic catalysis mechanism (see, e.g., 
Refs. [51, 52]). 4 



4 A very similar argument applies to a different observable as well. Recently it was shown that the QCD magnetization 
separates into spin- and orbital angular momentum-related contributions, and the spin term has been determined on the 
lattice [53]. In fact, the spin-related contribution can be written as the derivative of A/ vac ' r for a spin-1/2 quark, with respect 
to the gyromagnetic ratio g (replace a = 1/2 — s z with a — 1/2 — g/2 ■ s z in the calculation of Subsec. 3.3). The renormalization 
prescription of Eq. (3.15), on the other hand, contains the fixed <?* = 2 and, thus, does not contribute to the derivative with 
respect to g. Again, the result is of 0(B 2 ), in agreement with Ref. [53] - in contrast to the total vacuum free energy density 
that we calculate here, which is always of the form + 0(B 4 ). 
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There is a further consequence of the separation, Eq. (3.20). Let us consider the case of a charged pion, 
with m = m n , for which the sigma term is calculated using the Gell-Mann-Oakes-Renner relation, 



m /„. = m q iptpB=o -> = 2 — ^ — . (3.24) 



2 ^ - m 9 /| 

Inserting Eq. (3.21) into Eq. (3.20), we get, 

^ + 0{B A ) = A/ vac (m) + ^ Z g scalar (m). (3.25) 

Now we differentiate both sides with respect to the quark mass, and use the right hand side of Eq. (3.23). 
Expressing the condensate we obtain, 

= | . 3^) = 1 «^ 

2 o(m z ) om q 2 m f4 

where we inserted the expression (3.17) for the scalar QED renormalization constant, and employed Eq. (3.24). 
Using the value of ^scalar _ 1/(48^2^ resu lt reproduces the chiral perturbation theory formula [54] for 
the charged pion contribution to the condensate, up to 0(B 2 ). Altogether, this shows that the magnetic 
catalysis of the condensate is a direct consequence of the actual form of the electric charge renormalization 
constant in scalar QED. In particular, the fact that scalar QED is not asymptotically free (^calar > 0), 
implies that the quark condensate undergoes magnetic catalysis at T = 0, and increases (to leading order) 
quadratically with growing B. 

In the following, we will not consider the condensate, but only the free energy itself, for which it is not 
necessary to distinguish between m and m*. Moreover, we will exclude the pure magnetic term B 2 /2 from 
consideration, since it gives no information about the response of hadrons to the external field. Therefore, 
the vacuum pressure will be of 0(B 4 ). The superscript r denoting the renormalized free energy will also be 
dropped. 

3.5 Stability and spin channels 

We notice that the formula (3.19) is only well-defined for values of x = m 2 /2qB for which 

x+^-s>0 -> (2s - 1) qB < m 2 . (3.27) 

For any spin s > 1, this constraint gives a critical magnetic field B c , where the theory breaks down. Clearly, 
as B c is approached, the assumptions of the model - in particular, that hadrons are point-like particles - 
become incorrect. This implies that one is restricted to x > 1/2 for the s = 1 and x > 1 for the s = 3/2 
channel. The critical magnetic fields corresponding to the hadron is eB c (p ± ) = 0.6 GeV 2 . This 

observation forms the basis of the idea of a superconducting vacuum at high magnetic fields B > B c [55]. 

In Fig. 1 we plot the dependence (3.19) of the vacuum pressure Ap vac = — A/ vac on m 2 /2qB in the 
four spin channels s = . . . 3/2, with q/e = 1. We find that the vacuum pressure for the s = 3/2 channel is 
negative (and has been multiplied by —1 in the figure), whereas the other three are positive. In fact, above 
s = l, the sign of Ap v£LC starts to alternate, with spin-integer hadrons contributing positively and spin-half- 
integer hadrons negatively to the pressure. The s = 1/2 channel seems to be the only exception to this rule. 
The colored dots in Fig. 1 represent the corresponding lowest-lying particles at eB = 0.1 GeV 2 . Although 
for this magnetic field, the charged pion is seen to dominate, as the magnetic field grows, at some point the 
A 1 * 1 would clearly take over, turning the total vacuum free energy density / vac positive, and, accordingly, 
the T = pressure negative. This indicates an instability, which suggests that the HRG model in terms of 
the dispersion relation (3.2) is not applicable for s = 3/2. 
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Figure 1. Magnetic field dependence of the vacuum pressure in various spin channels for charge q/e — 1. The dots 
represent the corresponding values for a tt^ 1 (red), a proton (blue), a p ± (green) and a A* particle (gray) at a magnetic 
field eB — 0.1 GeV 2 . The dashed lines represent lower bounds on m 2 /2eB for each spin sector. The contribution 
from the A* has been multiplied by —1. 



4 Results 

We will make use of the thermodynamic relations (2.2)-(2.5) to determine the equation of state for nonzero 
magnetic fields. At B = 0, it is customary to normalize the free energy density, or, the pressure, by T 4 . As 
the magnetic field is switched on, the T = pressure is in general nonzero, and thus p/T 4 diverges as T 
decreases. We remark that the pressure at T = is nonzero even after dropping the pure magnetic energy 
B^./2 in Eq. (3.19), as (qB) 4 contributions are still present. Therefore, in the following we will plot the 
pressure and other EoS-related observables in physical units, without a normalization by T 4 . In Fig. 2, the 
individual contributions of various hadrons to the pressure are shown as functions of the temperature, for 
B = (left panel) and for eB = 0.2 GeV 2 (right panel). In the low temperature region (T < 100 MeV), 
the pressure is dominated by pions for vanishing magnetic field. As B increases, this dominance is lost, as 
vacuum contributions in the other spin channels - most importantly, for the hadron - arise. 




50 100 150 50 100 150 

T (MeV) T ( M eV) 

Figure 2. Individual contributions to the HRG pressure as a function of the temperature for eB = (left panel) 
and eB = 0.2 GeV 2 (right panel). Note the change in the contribution of charged particles, especially and p^, 
between the two panels. 
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Since, roughly speaking, the effective mass of charged particles is m e s = m 2 + qB{l — 2s) (see Eq. (3.2)), 
it increases with B for s = hadrons, but decreases for s = 1 particles. Accordingly, the thermal part of 
the pressure - which contains the Boltzmann weights exp(— m e ^/T) - is larger for p , whereas it smaller 
for ir^, as compared to the B = case. This effect is also visible in the temperature dependence of the tt^ 
and p^ contributions in the right panel of Fig. 2. We note moreover, that in our approach neutral particles 
are not affected by the magnetic field, since their gyromagnetic ratios are set to zero. 
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Figure 3. The equation of state in the HRG model. Shown are (from left to right and downwards) the pressure, the 
magnetization, the energy density, the entropy density and the speed of sound squared as functions of the temperature, 
for eB = (solid red lines), eB = 0.2 GeV 2 (dashed blue) and eB = 0.3 GeV 2 (dot-dashed green). 
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We proceed by performing the sum over hadrons to determine the dependence of thermodynamic ob- 
servables on B and T. These observables are the pressure p, the energy density e, the entropy density s, 
the magnetization ms and the speed of sound c s (for their definition, see Sec. 2). In Fig. 3, these quantities 
are plotted as functions of the temperature for eB = 0, 0.2 GeV 2 and 0.3 GeV 2 . The pressure is increased 
by the magnetic field in the hadronic sector, and, correspondingly, the magnetization is positive, indicat- 
ing a paramagnetic hadronic phase - with interesting implications regarding the deconfinement transition 
temperature, based on large N c arguments [56]. 

Since the entropy is written as the derivative of p with respect to T, it is insensitive to the vacuum 
contribution, and vanishes at zero temperature. Moreover, even at nonzero temperatures, s barely changes 
with B, in the range of magnetic fields under consideration. A much more pronounced signal can be seen 
in the speed of sound, which exhibits a dip, moving towards lower temperatures as B is increased. The 
minimum position of c 2 is one possible definition of the deconfinement transition temperature T c (see, e.g. 
Ref. [16]). The behavior of the speed of sound thus indicates that T c decreases as the external magnetic 
field grows, in agreement with the recent lattice results [12, 13]. 

5 Conclusions 

In this paper, we have developed a Hadron Resonance Gas model to study the QCD equation of state for 
nonzero magnetic fields. Using the renormalization properties of the free energy density, we have derived 
the relation Eq. (3.26), which connects the scalar QED /3-function and the mass-dependence of the quark 
condensate. This correspondence explains the well-known magnetic catalysis phenomenon, in terms of the 
renormalization group running of the scalar QED coupling, thereby relating two, seemingly very different 
concepts. 

We proceeded by investigating the individual contributions from hadrons to the pressure, and observed, 
that pions no longer dominate the low-temperature region if the magnetic field exceeds B > 0.2 GeV 2 , 
thus, it is essential to take into account higher-lying resonances - especially the hadron. By summing 
up the total pressure and considering its S-dependence, the magnetization is determined to be positive, 
showing that the QCD vacuum at low temperatures is paramagnetic. Using the magnetization and the 
pressure, the whole equation of state is reconstructed. The behavior of the speed of sound suggests, that the 
deconfinement transition temperature is lowered as the magnetic field is increased. We stress that the results 
obtained in the HRG model are reliable only at low temperatures and low magnetic fields. The limitation in 
the temperature is obviously given by the transition to the quark-gluon plasma phase, and as lattice results 
show, the HRG approximation at B = is reliable up to T ~ 150 MeV. Concerning the magnetic field, 
eB < m 2 ~ 0.6 GeV 2 must be fulfilled for the description to be consistent. The lattice determination of the 
EoS at nonzero magnetic fields is necessary to test the limitations of the HRG approach, and to confirm our 
findings. 
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A Formulae 



In this appendix we summarize the formulae that were used in the derivation of the vacuum free energy 
density in Subsec. 3.3. 

• For integration in d dimensions we use (see, e.g., Ref. [57]) 

r ^ v^tjp = 1 n-i/2-d/2) /2+d/2 

J^ 00 (2ny VP +M (4*)*/a r(-l/2) {M ) ' ^ 

and r(-l/2) = -2-v/tt. 

• The expansion of the T function around some negative integers is given by 

r(-l + e/2) = - 2 + 7 - 1 + O(e), F(-2 + e/2) = 1 - 1 + ] + O(e), (A.2) 

e e 2 4 

where 7 is the Euler constant. 

• The Hurwitz £ function is defined as 



k=0 

with the expansion [58] 



00 _. 

E^W = «*•*>• (A - 3) 



«-l + e/2, 1) « -i - ^ + I + |C(-l,x) + £>( e 2 ), (A.4) 



and asymptotic behavior [58], 



C'(-l, x) = 1 - ^ + f 1 - ^ + ^) log(x) + Oix-\ (A.5) 



12 4 V 12 2 2 

Writing the same integral in dimensional regularization with parameter e and scale fi, and in cutoff 
regularization with cutoff A, 



/ s e \ 



(A.6) 



we note that the two schemes are related as 



2 , / m 2 \ , / m 2 
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